Fragmentation of the photoabsorption strength in neutral and charged metal 

microclusters 



C. Yannouleaa3 and R.A. Broglia 

Dipartimento di Fisica, Universita di Milano, and INFN, Sez. Milano, 
1-20133 Milano, Italy, and The Niels Bohr Institute, DK-2100 Copenhagen 0, Denmark 

M. BradJ] 

The Niels Bohr Institute, DK-2100 Copenhagen 0, Denmark 
P.-F. Bortignon 

Istituto di Ingegneria Nucleare, CESNEF Politecnico di Milano, Italy, and INFN, LNL, Legnaro, Italy 
(Dated: 16 December 1988; Phys. Rev. Lett. 63, 255 (1989)) 

The line shape of the plasma resonance in both neutral and charged small sodium clusters is 
calculated. The overall properties of the multipeak structure observed in the photoabsorption cross 
section of spherical Nas and Na2o neutral clusters can be understood in terms of Landau damping. 
Quantal configurations are shown to play an important role. In the case of charged Na^J" and Na^ 
clusters a single peak is predicted that carries most of the oscillator strength. 
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Recent photoabsorption experiments [l[ have revealed 
the surface plasmon in small neutral sodium clusters. A 
qualitative overall account of the dependence of the res- 
onance frequencies on the cluster size can be achieved in 
terms of the extended ellipsoidal shell model 0, 0] and 
of the experimental static polarizabilities. The reported 
line broadening seems to obtain, as suggested in ref. [![, 
an important contribution from the coupling of the dipolc 
resonance to quadrupole shape fluctuations of the cluster 
0,1. 

Due to the fact that no complete photoabsorption 
curve has yet been experimentally determined, the de- 
tailed shapes and widths of the resonances are still rather 
uncertain. In particular, the contribution of direct plas- 
mon decay mechanisms to the line widths must be inves- 
tigated, especially in the case of small clusters, in view 
of the detailed microscopic calculations carried out by 
Ekardt [|| and Beck Q within the framework of the time- 
dependent local density approximation (TDLA). 

In the present paper, we study the structure of 
the plasmon resonance within the mean-field framework 
making use of the random phase approximation (RPA). 
Because these studies can most clearly be done in the 
case of spherical clusters, we presently restrict our analy- 
sis to the neutral Nas an d Na2o and to the charged Na^J" 
and Na^. It will be concluded that, in the case of the 
neutral Nas and Na2o, Landau damping is important and 
that the observed line shapes are the result of a detailed 
interplay between quantum size effects and residual in- 
teractions among particlc-hole excitations, as well as of 
the thermal fluctuations of the cluster shapes. In the case 
of the charged Nag and Na^ , Landau damping is prac- 
tically absent and the dipole strength essentially consists 
of a single peak exhausting most of the plasmon oscillator 
strength. 



The calculations to be discussed below parallel those 
carried out in studies of giant resonances in nuclei (cf., 
i.e. [8[). A discrete particle- hole basis of dimension Af is 
constructed and the Hamiltonian 



H = H a 



V, 



(1) 



sum of the Hatree-Fock Hamiltonian Ho and the residual 
interaction V, is diagonalized using the RPA. In what 
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FIG. 1: Self-consistent potentials [l(j associated with Na2o 
and Na^ and used in the solution of the single-particle Hamil- 
tonian Ho appearing in eq. {T}. For Na2o, the resulting single- 
particle levels relevant to the discussion of the dipole reso- 
nance, as well as some of the associated unperturbed particle- 
hole transitions, are also shown. 
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follows, we shall specify the single-particle potentials self- 
consistently in the spherical jellium-background model Q 
using the density variational formalism in a semiclassical 
approximation Q. These potentials are displayed in 
Fig. 1. 

The residual two-body interaction V is given by 



V(|ri-f 2 |) 



dV xc [p] 



Sdn-ril) . (2) 



1^1 — dp 
Here V^ c [p] = d£ xc [p]/dp is the exchange-correlation po- 



tential in the ground state. As in refs. we use the 

exchange-correlation energy density £;e C [p] of Gunnarsson 
and Lundqvist 

The RPA equations (see, e.g., ref. [13]) 



A B 

B* A* 



X n 



E,, 



-Yn 



are written in terms of the angular-momentum coupled 
matrix elements of the interaction ([2|), i.e. 



J 



where 



A(ph,p'h') - {e np . ip - e nhtlh )5 lptlp ,S lh!lh ,6 npinp ,6 nhtnh , = 2 R(ph',hp') (-) l » +l p' 

[(2l p + l)(2l h + l)(2l p , + l)(2l h , + lp ( l h A l p \ ( l h , X l pl 



(2A + 1) 











R A B(ph,p>h>) 



R(ph',hp')= / r 1 drir 2 dr 2 Tln p ,i p (ri)TZn h ,,i h Xr2)V(r 1 ,r2;X)U nh ,i h (r 1 )TZ np ,< ,(r 2 ) , 



r 



(3) 



(4) 



and where TZni,h(r) is the radial part of single-particle 
wave functions. The radial contribution of the two-body 
interaction ([2]) in multipole order A is given by 



V( ri ,r 2 ;X) = 



„X 



dV xc [p] 5(n-r 2 ) 2A + 1 



„A+i 



dp 



4tt 



where r< = min(ri,r2) and r> = max(ri, r 2 ). 

The indices appearing in eqs. © and ^ denote 
the number of nodes for the corresponding single-particle 
states. The orbital angular momenta of the particles and 
holes participating in the excitations are denoted by 
the total angular momentum of the excitation being A, 
which in the present calculation is set equal to 1 (dipole 
vibration). The 3j symbols appearing in take proper 
care of the angular momentum coupling, as well as of the 
parity conservation conditions. The factor 2 accounts for 
the spin degeneracy. 



The RPA eigenvectors are written as a linear combi- 
nation of particle-hole excitations in terms of the for- 
wardsgoing and backwardsgoing amplitudes, X n (ph) and 
Y n (ph) respectively, according to 

|n) = Y,[X n (j>h)\(ph- l ) Xil ) 

ph 

-{-) X+ "Y n {ph)\{hp^) x _J] , (5) 

where a singlet spin configuration is implied. 

The dipole transition probabilities associated with the 
state ([5]) can be written as 



B(El,Q^n) = - |H|A4(£1)||0>| 2 , 



where 



(n||M(£l)||0) - J2 [ x *n{ph; 1) + R A Y; L {ph- 1)] ( P \\M(El)\\h) 

ph 



(6) 



r 



are reduced matrix elements [13| of the dipole operator 
M{El;p) = ^j3ery lfl (f). 



The radial wave functions R ni ,ii are calculated by di- 
agonalizing the single-particle Hamiltonian Hq in a ba- 
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sis including AT = 25 harmonic oscillator major shells. 
(For the oscillator parameter of this basis, we have used 
fiLj = 1.1 eV. The results do not, however, depend on 
this choice.) 

The associated unperturbed particle- hole excitations, 
examples of which are displayed in Fig. 1, exhaust 
the Thomas-Reiche-Kuhn sum rule S(E1) = Nfi 2 e 2 /2 m, 
that is 



1.0 



ph 



, p - e nh , lh )\^2/3 (p\\M(El)\\h)\ 2 = S(E1) 



where e. ni ,ii are the single-particle energies. This result 
is also valid for the correlated eigenstates (© and associ- 
ated eigenvalues E n , since the RPA preserves the energy- 
weighted sum rule (EWSR) . 

In Fig. 2 we display the oscillator strength functions 
versus excitation energy for the RPA dipole in the case 
of the neutral Nas and Na2o clusters, as well as in the 
case of the charged Na^ and Na^ clusters. The unper- 
turbed oscillator strength for the neutral clusters is also 
displayed. 

A conspicuous feature of the RPA results for the neu- 
tral clusters is the sizeable amount of Landau damping. 
In fact, the variances of these response functions are 
a w 0.57 eV in both cases implying a ratio a/E « 0.25, 
E here being the energy centroids of the resonances. 
Nonetheless, the identification of the collective states is 
quite unique. Indeed, in the case of Nas, there is a single 
state at rj 2.8 eV exhausting w 75% of the EWSR, while 
about the same strength is distributed among two lines 
located at 2.6 eV and 2.9 eV in the case of Na2o- As 
befitted collective states, the associated wave functions 
exhibit several (about ten) forwardsgoing amplitudes (X- 
components) that are larger than 0.1. Also, most of the 
backwardsgoing amplitudes (Y-components) contribute 
constructively to the transition amplitude ([6]) . The wave 
functions associated with the rather weak states which 
are strongly red shifted with respect to the collective 
states are, on the other hand, dominated by a couple 
of components. We note that the results shown in the 
middle of Fig. 2 for Na2o display an extent of fragmen- 
tation similar to the TDLDA calculation by Ekardt (cf. 
Fig.6ofref. 0). 

The plasmon strength function is drastically modified 
in the case of the charged Nag and Na^ clusters, where 
a single peak lying at w 3.1 eV (Nag") and at f» 3.0 
eV (Na+) exhausts w 93% and w 83% of the EWSR, 
respectively. This result seems to be consistent with re- 
cent experimental findings reported [3] for the potassium 
clusters Kg and K^. 

The difference between positively charged and neutral 
clusters is due to the difference in the associated aver- 
age potentials (cf. Fig. 1). Indeed, the potentials for 
the charged clusters are deeper than the potentials for 
the corresponding neutral clusters (observe that the po- 
tential for the neutral clusters exhibit an almost constant 
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Oscillator strength function for the photoabsorption 



Na 20 , Na;}", 



and Najj 



clusters. In the upper figures, 



the unperturbed oscillator strengths for the neutral clusters 
is displayed. Also, in the upper figures, the strength values 
above 3 eV have been multiplied by 100. 



central depth at w —5.4 eV). For Nag" the increase of the 
potential depth at the center amounts to ~ 56%, while 
for Na^ the corresponding increase is w 42%. 

It is interesting to compare the moments of the dipole 
strength function obtained in the present microscopic 
RPA calculation for the case of neutral clusters with 
those obtained semiclassically in ref. [l(J. The cubic 
energy-weighted sum rules and the centroids E agree 
within 2 — 4%; the static dipole polarizabilites a (given 
by the negative energy-weighted sum rules) agree within 
4 — 10%. This offers support for the ability of the semi- 
classical method proposed in ref. to predict the av- 
erage properties of the mean-field response function. 

In order to compare the present calculations with the 
experimental data, one has to take into account the cou- 
pling of the electronic dipole oscillations to thermal sur- 
face fluctuations of the whole cluster 0, [|| . This is done 
only for the case of neutral clusters for which the most 
systematic body of information is available at present [l5j | . 
In keeping with the findings of ref. [Bj], this can be ap- 
proximately simulated by folding the peaks of the RPA 
spectrum with Lorentzian functions displaying a damp- 
ing factor T/huJo « 0.1 and 0.06 for Nag and Na2o, re- 
spectively. The results are shown in Fig. 3 in comparison 
with the experimental data [l6j |. Aside from a shift of 
« 10%, the overall structure observed in these two cross 
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FIG. 3: Photoabsorption cross sections per atom (for Nas 
and Na2o) resulting from the folding of the RPA oscillator 
strengths with Lorentzian shapes normalized to unity are 
shown (solid line). The widths of the Lorentzians are speci- 
fied by T/huio ~ 0.1 and 0.06 for Nag and Na2o, respectively 
[jj, fiu)o being the peak energy of a given Lorentzian. The 
experimental data[T(| are also shown. The dashed line is the 
solid line shifted so that there is a maximum overlap with the 
data. 



sections is well reproduced. In particular for Na2o, the 
dip seen in the experimental cross section near 460 nm, 
as well as the weak bump in Nas near 590 nm, seem to 
come out of the present RPA description. 

The shift of the calculated dipole peaks with respect to 
the observed values for the neutral clusters is consistent 
with the fact that the predicted static dipole polarizabil- 
ities within the LDA model [UHU 12 1 are lower than the 
experimental ones by sa 20%. This is because the model, 
when applied to the exchange part of the Coulomb inter- 
action between the electrons, provides too much screen- 
ing, and thus fails to produce for neutral clusters the 
expected asymptotic 1 jr Coulomb behavior of the mean 
field. As a consequence, the tail of the electron densities 



(the so-called 'spill-out') is underestimated; since the lat- 
ter is known to be correlated to both the dipole polariz- 
abilities 17 1 and the surface plasmons [1(3], these quan- 
tities are also underestimated. Indeed, a 'self-interaction 
correction' aimed at a better treatment of the Coulomb 
exchange within the Kohn-Sham approach [l8j can im- 
prove the values for the dipole polarizabilities. It remains 
an open question to which extent an improved treatment 
of the Coulomb exchange will affect the positions of the 
unperturbed particle-hole excitations, and thus the de- 
tailed structures of the photoabsorption cross sections 
obtained in the present investigation. 

It can be concluded that direct decay of the giant 
dipole resonance in small neutral sodium clusters plays 
an important role in the fragmentation of the associated 
strength function, a process that seems to be essentially 
absent in the case of the positively charged clusters. The 
actual location and potential fragmentation of the dipole 
peak in Nag and Na2o is the result of a delicate balance 
between shell structure and residual interaction, which 
might be further unravelled by carrying photoabsorption 
measurements at low temperatures. A detailed mapping 
of the strongly red shifted dipole peaks and of their struc- 
ture could provide some crucial tests of the exchange 
(plus correlation) part of the effective interaction between 
the electrons in metal clusters. 
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